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We present simple models of particulate materials whose mechanical in- 
tegrity arises from a jamming process. We argue that such media are 
generically "fragile", that is, they are unable to support certain types 
of incremental loading without plastic rearrangement. In such models, 
fragility is naturally linked to the marginal stability of force chain net- 
works (granular skeletons) within the material. Fragile matter exhibits 
novel mechanical responses that may be relevant to both jammed colloids 
and cohesionless assemblies of poured, rigid grains. 



1 Introduction 

In this paper, we consider the relation between jamming (a kinetic process) 
and the laws of static mechanical equilibrium of particulate media. We describe 
first a simple model of jamming in a colloid, sheared between parallel plates. 
We assume that force chains (linear arrays of rigid particles in contact) develop 
along the major stress compression axis and span the sample. The resulting 
jammed state can support a shear stress indefinitely, but does this by virtue of 
a self-organized internal structure (the force chain array) which has developed 
in direct response to the applied load itself. If a different load is now applied 
(e.g. if the material is sheared in some different direction), the force chains 
cannot sustain the load but must flow and rejam in a different configuration. 
The model thus provides a simple example of a "fragile" material - one which 
cannot support certain types of infinitesimal stress increment without plastic 
reorganisation. 
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Fig. 1. (a) A jammed colloid (schematic). Black: force chains; grey: other 
force-bearing particles; white: spectators, (b) Idealized rectangular network of force 
chains. 

The continuum mechanics of fragile materials is very different from conven- 
tional elastic or elastoplastic theory. We propose that fragility may be generic 
in jammed systems, and argue further that a pile of cohesionless poured sand 
can be viewed as jammed in the required sense: its mechanical integrity results 
purely from the applied load (gravity). Some simple models that we have re- 
cently developed for poured sand such as the fixed principal axes (fpa) 
model, indeed display fragile behaviour; in fact they correspond to assuming 
a particular simplified geometry for the granular skeleton of force chains. 

Fragile models require that the behaviour of a pile of cohesionless grains (frag- 
ile) is quite different from a hypothetical pile where each grain is firmly glued 
to its neighbours on first coming to rest (elastic). This is plausible because for 
random packings of glued grains, a finite fraction of the interparticle forces 
will be under tension (forbidden in the cohesionless case) and will remain so 
until the load is large enough to cause appreciable elastic deformation of indi- 
vidual grains. For simple fragile models such as FPA, the same criterion defines 
a crossover between fragile behaviour and a form of (anisotropic) elastoplas- 
ticity. A fuller discussion of these ideas is given in Q, where further references 
may be found. 



2 Jamming in Colloids: A Simple Model 

Consider a concentrated colloidal suspension of hard particles, confined be- 
tween parallel plates at fixed separation, to which a shear stress is applied. 
Above a certain threshold of stress, this system exhibits enters a regime of 
strong shear thickening; see, e.g. ||. The effect can be observed in the kitchen, 
by stirring a concentrated suspension of corn-starch with a spoon. In fact, com- 
puter simulations suggest that, at least under certain idealized conditions, the 
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Fig. 2. (a) A force chain of hard particles (any shape) can statically support only 
longitudinal compression, (b) Finite deformability allows small transverse loads to 
arise. 

material will jam completely and cease to flow, no matter how long the stress 
is maintained ||. In these simulations, jamming apparently occurs because the 
particles form "force chains" along the compressional direction (Fig. p] (a)). 
Even for spherical particles the lubrication films cannot prevent direct inter- 
particle contacts; once these arise, an array or network of force chains can 
indeed support the shear stress indefinitely. (We ignore Brownian motion, 
here and below, as do the simulations; this could cause the jammed state to 
have finite lifetime.) 

To model the jammed state, we start from a simple idealization of a force chain: 
a linear string of at least three rigid particles in point contact. Crucially, this 
chain can only support loads along its own axis (Fig.0 (a)): successive contacts 
must be collinear, with the forces along the line of contacts, to prevent torques 
on particles within the chain ||. Note that neither friction at the contacts, 
nor particle aspherity, can change this "longitudinal force" rule. 

As a minimal model of the jammed colloid, we take an assembly of such 
force chains, characterized by a unique director n, in a sea of "spectator" 
particles, and incompressible solvent. This is obviously oversimplified, for we 
ignore completely any collisions between chains, the deflections caused by 
weak interactions with the spectator particles, and the fact that there must 
be some spread in the orientation of the chains themselves |J. With these 
assumptions, in static equilibrium, with no body forces acting, the pressure 
tensor (defined as Pij — a^, with the usual stress tensor) must obey 



Here P is an isotropic fluid pressure, and A (> 0) a compressive stress carried 
by the force chains. 

Eq. (Q) defines a material that is mechanically very unusual. It permits static 
equilibrium only so long as the applied compression is along n; while this 



Pij = Pdij + A UiHj 
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remains true, incremental loads (increase or decrease in stresses at fixed major 
compression of the stress tensor) can be accommodated reversibly, by what 
is (at the particle contact scale) an elastic mechanism. But the material is 
certainly not an ordinary elastic body, for if instead one tries to shear the 
sample in a slightly different direction (causing a rotation of the principal stress 
axes) static equilibrium cannot be maintained without changing the director 
n. Now, n describes the orientation of a set of force chains that pick their 
ways through a dense sea of spectator particles. Accordingly n cannot simply 
rotate; instead, the existing force chains must be abandoned and new ones 
created with a slightly different orientation. This entails dissipative, plastic, 
reorganization, as the particles start to move but then re-jam in a configuration 
that can support the new load. 



3 Jamming and Fragile Matter 



Our model jammed colloid is an idealized example of "fragile matter": it can 
statically support applied shear stresses (within some range), but only by 
virtue of a self-organized internal structure, whose mechanical properties have 
evolved directly in response to the load itself. Its incremental response can 
be elastic only to compatible loads; incompatible loads (in this case, those of 
a different compression axis), even if small, will cause finite, plastic reorga- 
nizations. The inability to elastically support some infinitesimal loads is our 
proposed definition of the term "fragile" (which, up to now, has not been given 
a precise technical meaning in this context). 

We argue that jamming may lead generically to mechanical fragility, at least 
in systems with overdamped internal dynamics. Such a system is likely to 
arrests as soon as it can support the external load; since the load is only 
just supported, one expects the state to be only marginally stable. Any in- 
compatible perturbations then force rearrangement; this will leave the system 
in a newly jammed but (by the same argument) equally fragile state. This 
scenario is related, but not identical, to several other ideas in the literature 
fl8Hl3f . These include the emergence of rigidity by successive buckling of force 



chains in glasses and granular matter ||; the concepts of self-organized crit- 
icality (soc) JlO| (see also ||), and those of mechanical percolation which 



underly recent "hypoplastic" models of granular matter [ H] . Fragility is also 
connected with recent ideas concerning isostaticity and marginal mechanics in 
(frictionless) sphere packings [0,0 (see Section 0). 

Our ideas are, in addition, reminiscent of the (much older) concept of a self- 
selecting critical state in theories of soil mechanics Jl4|]. However the latter is 
primarily concerned with dilatancy: the tendency of dense particulate media 
to expand upon shearing. Jamming can be viewed as the constant-volume 
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counterpart of this process: if expansion is prevented, jamming results. 



4 Two Types of Fragility 



Consider the idealized jammed colloid of (Fig. [I] (a)). So far we allowed for 
an external stress field (imposed a the plates) but no body forces. What body 
forces can it now support without plastic rotation of the director? Various 
models are possible. One is to assume that Eq. (U) continues to apply, with 
P(r) and A(r) now varying in space. If P is a simple fluid pressure, a localized 
body force can be supported only if it acts along n. Thus (as in a bulk fluid) no 
static Green function exists for a general body force. (Note that, since Eq. ([!]) 
is already written as a continuum equation, such a Green function would 
describes the response to a load that is localized in space but nonetheless acts 
on many particles in some mesoscopic neighbourhood.) 

For example, if the particles in Fig. [I] (a) were to become subject to a grav- 
itational force along y, then the existing force chains could not sustain this 
but would reorganize. Applying the longitudinal force rule, the new shape is 



easily found to be a catenary, as realized by Hooke ||15|| , and emphasized by 
Edwards ||. On the other hand, a general body force can be supported, in 
three dimensions, if there are several different orientations of force chain, pos- 
sibly forming a network or "granular skeleton" [{7,11,17,18 . A minimal model 
for this is: 

Pij = Ai riiUj + A 2 mirrij + A 3 klj (2) 



with n, m, 1 directors along three nonparallel populations of force chains; the 
A's are compressive pressures acting along these. Body forces cause A 12j 3 to 
vary in space. 

We can thus distinguish two levels of fragility, according to whether incompat- 
ible loads include localized body forces (bulk fragility, e.g. Eq. (P), or are lim- 
ited to forces acting at the boundary ( boundary fragility, e.g. Eq. (^)). In dis- 
ordered systems one should also distinguish between macro-fragile responses 
involving changes in the mean orientation of force chains, and the micro-fragile 
responses of individual contacts. Below we focus on macro-fragility, but in 
practice the distinction may become blurred. In any case, these various types 
of fragility should not be associated too strongly with minimal models such 
as Eqs. (110). It is clear that many granular skeletons have a complex network 
structure where many more than three directions of force chains exist. Such a 
network may nonetheless be fragile; see Section [7] below. 



5 



5 Fixed Principal Axis (fpa) Model 



Returning to the simple model of Eq. (|), the chosen values of the three 
directors (two in 2-d) clearly should depend on how the system came to be 
jammed (its "construction history"). If it jammed in response to a constant 
external stress, switched on suddenly at some earlier time, one can argue 
that the history is specified purely by the stress tensor itself. In this case, if 
one director points along the major compression axis then by symmetry any 
others should lie at rightangles to it (Fig. [I] (b)). Applying a similar argument 
to the intermediate axis leads to the ansatz that all three directors lie along 
principal stress axes; this is perhaps the simplest model in 3-d. One version 



of this argument links force chains with the fabric tensor [|TTJ, which is then 
taken coaxial with the stress 



With the ansatz of perpendicular directors as just described, Eq. (g) becomes a 
"fixed principle axes" (fpa) model 0,0]. Although grossly oversimplified, this 
leads to nontrivial predictions for the jammed state in the colloidal problem, 
such as a constant ratio of the shear and normal stresses when these are varied 
in the jammed regime. Such constancy is indeed reported by Laun || in "the 
regime of strong shear thickening"; see H. 



6 Granular Materials 



We now turn from colloids to granular materials. Although the formation of 
dry granular aggregates under gravity is not normally described in terms of 
jamming, it is a closely related process. Indeed, the filling of silos and the 
motion of a piston in a cylinder of grains both exhibit jamming and stick-slip 
phenomena associated with force chains; see [|19|]. And, just as in a jammed 
colloid, the mechanical integrity of a sandpile disappears as soon as the load 
(in this case gravity) is removed. 

In the granular context, a model like Eq. (|2]) is interpreted by saying that a 
fragile granular skeleton of force chains is laid down at the time when particles 
are first buried at the free surface; so long as subsequent loadings are com- 
patible with this structure, the skeleton will remain intact. If in addition the 
skeleton is rectilinear (perpendicular directors) this forces the principal axes 
to maintain forever the orientation they had close to the free surface (fpa 
model). However, we do not insist on this last property and other models, 
based on an oblique family of directors, have also been developed |]|f]. (The 
construction history of a sandpile allows this since the orientation of the free 
surface and/or gravity provides a reference direction in addition to that given 
by the stress tensor at the onset of jamming [|TJ.) Note also that, for a conical 
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sandpile we require in addition to Eq. (0) one further relation among stresses, 
found for example by assuming that A 2 = A 3 everywhere; see [0. It turns out 
that the fpa models account quite well for the forces measured experimentally 
beneath conical piles of sand, constructed by pouring cohesionless grains from 
a point source onto a rough rigid support [|20|JT|] . (Note that the two dimen- 
sional wedge - may have special features and the effectiveness of fpa 
for this is much less clear ||.) 

As mentioned previously, fragile models such as this show very different me- 
chanics from conventional forms of elasticity or elastoplasticity. For example, 
in 2-d, when combined with stress continuity (diO~ij = pgj for sand under 
gravity), Eq. (Q) gives differential equations for the stress tensor which are 
hyperbolic 0. With a zero-force boundary condition at the upper surface of 
a pile, this gives a well-posed problem: the forces acting at the base follow 
uniquely from the body forces by integration. (Analogous remarks apply to 
Eq. (0) in 3-d.) If different forces are now imposed at the base, rearrangement 
is inevitable ||21|| . (This is boundary- fragile behaviour.) The same does not hold 
j|] within a traditional elastoplastic modelling approach [^2j whose equations 



are elliptic in elastic zones. In such models of the sandpile the forces acting 
at the base cannot be found without specifying a displacement field there. 
To define this displacement, one would normally invoke as reference state the 
one in which the load (gravity) is removed. For cohesionless poured sand, this 
state is undefined [23|, just as it is for a jammed colloid which, in the absence 



of the applied shear stress, is simply a fluid. 

One route to an elastic reference state is to consider a hypothetical sandpile 
where each grain becomes firmly "glued" to its neighbours (or the base), upon 
first coming to rest. The resulting medium is surely elastic, and must therefore 
be governed by elliptic equations. This does not mean that it is a conventional 
homogeneous elastic continuum (for which the states of zero strain and zero 
stress coincide). Indeed, a glued pile built under gravity will certainly have 
nonzero stresses if gravity is later removed HJ. More importantly, for a typical 
disordered packing of near-rigid, glued grains, there will arise many tensile 
contacts even under a purely compressive external load. Thus the problem of 
glued and unglued piles might, in practice, be extremely dissimilar. 



7 Isostaticity, Marginal Packings, and Fragility 



The problem of cohesionless granular media is a highly nonlinear one, because 
(a) grains are typically very rigid (on a scale set by the stresses that arise), 
and (b) there can be no tensile forces between any grains in the entire system. 
The first issue is addressed in |4| where we show that, for some fragile models, 
there is in fact a smooth crossover to more conventional elastoplastic physics 
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when grains are deformable, with fragility emerging as the limiting behaviour 
for rigid particles. 



In this rigid particle limit, where the longitudinal force rule (lfr) becomes 
valid, fragility will be recovered in a granular skeleton of force chains whose 
coordination number z = 2d with d dimension of space (e.g. Fig. p] (b) in two 
dimensions). This is the same rule as applies for packings of frictionless hard 
spheres, which also obey the lfr - not because of force chains, but because 
there is no friction. Indeed, regular packings of frictionless spheres, which 
show "marginal" or "isostatic" mechanics have been studied in detail recently 



12| . Note that although our idea of fragility is closely related to the concept of 
isostaticity (see e.g. |U|), it is apparently not identical; in our models of fragile 
materials the isostatic condition applies only to a loadbearing substructure 
(the force skeleton) and not the whole packing. It would be interesting to test 
this idea numerically by measuring the mean coordination number among the 
subset of particles that carry strong forces. 

Interestingly, Moukarzel has recently argued jT3] that random packings of 
frictionless spheres generically become isostatic [z = 2d) in the rigid parti- 
cle limit. His arguments appear to depend only on the lfr and the absence 
of tensile forces, so they should, if correct, equally apply to any granular 
skeleton that is made of force chains of rigid particles. It remains to be seen 
whether Moukarzel's (somewhat unintuitive) arguments can be made rigorous. 
Pending this, we prefer for the moment our own qualitative reasoning, which 
generically links fragility to jamming, as the motivation for setting z = 2d 
in our simplified models (such as fpa) of the granular skeleton. Such models, 
in which the skeleton is represented as a rectilinear or oblique array with no 
disorder, are obviously convenient for calculation. However, disorder will not 
remove the fragility (though it will change Eqs. fll]]^)) unless it causes the 
mean coordination number of the skeleton to increase. 



8 Conclusions 



The jammed state of colloids, if it indeed exists in the laboratory, has not 
yet been fully elucidated by experiment. It is interesting the even very simple 
models such as Eq. ([!]) can lead to nontrivial and testable predictions (such 
as the constancy of certain measured stress ratios). Such models suggest an 
appealing conceptual link between jamming, force chains, and fragile matter 
H. However, further experiments are needed to establish the degree to which 
they are useful in describing real colloids. 

For granular media, the existence of tenuous force-chain skeletons is clear 
|7| ]TI| , |IE| -|IB| ; the question is whether such skeletons are fragile. Several theo- 
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retical arguments have been given, above and elsewhere, to suggest that this 
may be the case, at least in the limit of rigid particles. Moreover, simula- 
tions show strong rearrangement under small changes of compression axis; the 
skeleton is indeed "self-organized" [fL7|J18|l . Experiments also suggest cascades 



of rearrangement [0,|R]] in response to small disturbances. These findings are 
consistent with the fragile picture. 

The mechanics of fragile models such as Eqs. ([i]j2|) differ strongly from those of 
conventional elasticity. For example, if an infinitesimal point force is exerted 
on top of a layer of sand, we expect the resulting pressure distribution at the 
base to form an annulus 0; an elastic model would predict the maximum 
pressure to be directly beneath the applied force. It could be very fruitful to 
perform new (but quite simple) experiments of this kind. 
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